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ABSTRACT 

CNI \ We determine a crucial feature of the dark halo density distribution from 

the fact that the luminous matter dominates the gravitational potential at 

> 

■ about one disk scale-length i?£>, while at the optical edge R op t — 3i?d the 

OO ' 
CO 

the kinematics of 137 spirals we find that the DM halo density profiles are 
. self-similar at least out to R op t and show core radii much larger than the cor- 

o ■ 

■ responding disk scale- lengths. The luminous regions of spirals consist of stellar 

■ 

Oh' disks embedded in dark halos with roughly constant density. This invariant 

DM profile is very difficult to reconcile with the fundamental properties of 
c/3 , the density distribution of CDM halos. With respect to previous work, the 

present evidence is obtained by means of a robust method and for a large and 



complete sample of normal spirals. 

1 INTRODUCTION 

Rotation curves (hereafter RCs) of spiral galaxies do not show any Keplerian fall-off: this 
implies the presence of an invisible mass component (Rubin et al. 1980; Bosma 1981; see 
Salucci and Persic, 1997). More precisely, the mass distribution of stars and gas does not 
match that of the gravitating matter (Persic and Salucci 1988, Persic, Salucci and Stel, 1996, 
hereafter PSS, see also Corbelli and Salucci, 1999); the discrepancy increases with increasing 
radius and, at a given radius measured in units of disk length-scales Rd, it increases with 
decreasing galaxy luminosity (Persic & Salucci 1988, 1990a,b; Broeils 1992). 

Each individual circular velocity V(R) and the spiral Universal Rotation Curve can both 
be represented, in the optical regions, by a linear law: (Rubin et al. 1980, see Persic & Salucci 
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(1991) and PSS for details) 



V(R) ~ V opt [l + V (-^- - l)] 0.4< R/R opt < 1.2 (la) 



where -R OJrf = 3.2i? d , V opt = V(R opt ), and 



V(^) = 0.10-1.35(/o 9s ^) (11) 



In (lb) the range is -0.1 < V = dlogy(R) 



cflogi? 



< 0.7 and the r.m.s. is 0.05 (see PSS and Fig 

Ropt 



!))• □ Notice that we will freely interchange luminosities and V opt = V(R opt ), given their 
high degree of correlation. 

The analysis of the URC and/or of individual RC has provided crucial knoweledge on the 
main global properties of the dark and luminous matter (e. g. Salucci & Persic, 1997). On the 
other hand, to invesatigate the local properties of the dark halos (e.g. the central density) 
is quite difficult, especially for luminous matter (LM) dominated objects (van Albada et al 
1985, but see the exception of high-resolution RC's, Borriello and Salucci, 2000). However, 
as a result of recent substantial observational and theoretical progress, a proper investigation 
of the spiral's halos mass structure is now possible. In fact, from the study of a large number 
of high-quality RC's recently available (Persic & Salucci, 1995) it has emerged that: 
i) the dark matter follows a regime of Inner Baryon Dominance (IBD) according to which 
in every normal spiral there is a transition radius Ribd 



inside which the luminous matter totally accounts for the whole mass distribution (see 
Ratnam and Salucci 2000, Salucci and Persic 1999ab, Salucci et al, 2000). This allows us 
to address the issue of the degeneracy problem raised by Van Albada, et al. (1985): the 
URC (and individual rotation curves) do show, in their profile, the kinematical signature of 
a transition between an inner LM-dominated region and an outer DM-dominated one (e.g. 
Salucci and Persic, 1999b). 

ii) dark halos are distributed very differently from the various "luminous" components (Cor- 
belli and Salucci, 1999). 

These findings allow us 1) to improve the determination of the disk mass to the level 
required for investigating the DM distribution and 2) to relate the dark halos around galaxies 
to collision-less non-baryonic cosmological structures. 

* We take R op t as the reference disk scale to follow PSS. We can use any other multiple of R d to specify the URC: no result 
of this paper changes. The URC for the whole range of available data, 0.1-R O pt < R < 2R op t, is given in PSS 
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The aim of this letter is to derive, for a large and complete sample of spirals, the density 
profile of the dark halos at the edge of the disks which are embedded in them and reveal dark 
constant-density regions of obvious cosmological importance. The evidence for core radii of 
DM halos, obtained so far by means of mass modeling of a (small) number of DM-dominated 
RC's (e.g. Flores and Primack 1994, Moore 1994, Burkert 1995), is being questioned in the 
light of the intrinsic uncertainty of the analysis itself (Burkert and Silk, 1997, van den Bosch 
1999). In many cases, in fact, the standard RC fitting method has difficulty in discriminating 
a NFW (Navarro, Frenk and White, 1996) density profile, which has Vh(R) oc R 1 ^ 2 in the 
center, from a const ant- density one, with Vh(R) oc R. 

We tackle this issue by resorting to the method of Persic and Salucci (1990b) in which 
the power law slope of the dark halo velocity at the disk edge is derived by means of a robust 
and straightforward procedure, which ultimately exploits the fact that, at ~ R op t, the dark 
halo is always the main density component, even when it is a negligible mass component at 
about Rd- \[ 

In detail, the goal of this letter is to detect a clear and reliable feature of the dark 
matter distribution, relevant on its own and (probably) at variance with the structural 
properties of standard CDM halos. Let us stress that the issue of establishing the actual 
CDM halo properties, or that of investigating whether non-standard CDM scenarios may be 
in agreement with observations, are beyond the scope of this work. 

The plan of this letter is the following: in section 2 we describe the RC sample and we 
derive the local DM slope Vh for all of the objects, in section 3 we compare the observed 
halo mass distributions with the CDM prediction and comment on our main results. 



2 THE HALO MASS DISTRIBUTION AT R OPT 
Let us start from the condition of rotational equilibrium: 

V\R) = V* m (R) + V*{R), (3a) 



where V^ m (i2) = Vg(R) + V^(R) + V 2 (R), with obvious notation, is the quadratic sum 



V (R) 

of the three "luminous" components: gas, disk, and bulge. We define (3 = yf^k 

V 2/d\ 



and 

Ropt 



V 2 (R) 



. In late type spirals, the exponential thin disk (of mass Ma) is by far the 

Ropt 



t This can be seen by combining eq (1) with the Poisson equation (see, e.g. Fall and Efstasthiou, 1980). 
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Figure 1. Rotation curve slope at R op t as a function of V op t- Triangles refer to dwarf spirals. 

main contributor to Vi um (R), at R ~ R opt (e.g. Verheijen, 1997; Rhee, 1997; Persic, Salucci 
& Ashman, 1993), and so 

P - Pd (36) 

with : 

_ _ (T n K„ - T, KA\, ... 

(4) 



Vl(x)=P d V^ 2{l0K0 ~ IlKl)W6X 



where x = R/R opt , I n ,K n are the modified Bessel functions of nth-order, and using the 
previous definitions, V^(R opt ), the DM velocity contribution at R opt is given by 

V2(R opt ) = (l-f3)V£ t (5) 

Let us define: 

^ d\ogV d {R) 



dlogR 



Ro 



V lurn, — 



d\ogVi um (R) 



dlogR 



Ropt 



(6a) 



(66) 



and V/i, V 9 , Vb in the same way of the l.h.s of equation (6a) when V d is substituted by 
Vh, V g and V& (the subscripts h,g,b refer to halo, gas and bulge respectively). In view of 
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the argument ii) in the previous section, we can identify Vh(R) with the contribution of a 
non-baryonic dark component. 
From the previous eqs. 

V,„ m = ( 1 ^r) V, + (^J^)\ + (t^t)^ (6c) 

v lum\ Jri 'opt) v lum\ ri 'opt ) v lum\ Jri 'opt) 

For a bulge-less gas-free spiral, 

V lum = V d = -0.273 (7a) 

from eqs. (4) and (6a). In this case, given its definition and the self-similarity of spiral 
disks, Vi um is strictly a constant. The contribution of bulge to V/ um can be totally neglected 
(see Persic, Salucci and Ashman, (1993) for details) while that of the gaseus disk can be 
evaluated in spirals with HI measurements: in a large sample (Rhee, 1997) we find that 
typically ( J^TK ) 2 ~3-6x 10" 2 , and V„ ~ 0.5 that leads to V imn ~ 0.9V d . For a 
sub-sample of the present sample we have (PSS): 

V lum = -0.24 ±0.03(2(7) (76) 

It is worth to stress that no result of this paper changes l)by assuming V; Mm according to eq 
(7a) rather than to eq (7b), or 2)by neglecting the (small) variance of Vi um among spirals. 
By differentiating eq. (3a) we arrive to (Persic & Salucci 1990b): 

which expresses the DM halo velocity slope in terms of the RC slope and of the LM mass 
fraction at Ropt- 

Note that, for high-quality RC's (SV < 0.05), the estimate of from eq (8) is very 
robust. In fact, in DM dominated objects, ((3 < 0.5), the uncertainties on f3 do not affect the 
estimate of Vh, while in LM-dominated objects, {(5 > 0.5), a reasonably good knowledge of 
f3: 5(3 < 0.1 suffices to estimate within a reasonable uncertainty (i.e. <5V/i < 0.2). 

Once we include the radial dependences of all the quantities, eq (8) is in principle valid 
at any radius. However, for R < R opt , (3{R) — > 1 very quickly, so that even a small error 
in (3 or V will strongly affect the estimate of Vh(-R). On the other hand, for R > R opt , the 
number of available RC's rapidly decreases. 

The sample of RC's we use in this letter is presented in PSS and Salucci and Persic (1997). 
It includes 131 rotation curves of spirals and 6 for dwarfs, all with a reliable profile out to 
R opt . This has been ensured by the following selection criteria: each RC (a) extends out to a 
R > Ropt] (b) has at least 30 velocity measurements distributed homogeneously with radius 
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Figure 2. The f3 vs V relationship for the present sample, (small filled circles) with their best fit (solid line). For a comparison 
we plot also the values derived from the URC and binned over V op t as in PSS (large filled cirles). 

and between the two arms; and (c) shows no asymmetries or non-circular motions. For 21-cm 
RCs we require that the beam-size should be < 1/27?^ in order to limit the uncertainties 
due to beam smearing. Incidentally, no LSB RC is found to satisfy these conditions, and 
this points to an intrinsic difficulty in studying these objects. 

The values of the V's are estimated by fitting with eq (1) each of the 137 RCs of the 
sample: the related uncertainty is small: <5V ~ 0.02 — 0.05; the values of V opt , V and <5V, 
are given in Tables 1 and 3 of PSS and shown in Fig (1) as function of V op t- It is worth to 
point out that a) the range of V along the luminosity sequence is large, ~ 0.7, while, the 
corresponding r.m.s., at a given luminosity is small < 0.1: the DM properties are likely to 
vary with galaxy global properties rather than random, b) we do not assume the URC (i.e. 
eq. (lb)): the values of V are derived from each RC . 

The LM fraction f3 is derived as in Salucci and Persic (1999a), Salucci et al. (1999) 
and Ratnam & Salucci (2000), i.e. by fitting the inner parts of each rotation curve, ( < 
R < Rjbd)i with (only) the circular velocity of an exponential thin disk given by Eq. 
(4). 0. In this region, the only-disk mass model reproduces, with no free parameters, the 
normalized rotation curve V(R)/V opt and, with a suitable choice of the parameter /3, the 
full curve V(R). This is shown for coadded RCs in Figure (3) and for individual RCs in 
Salucci et al., (2000) and in Ratnam and Salucci (2000)). The excellent match we always 

' In very few cases a bulge component has been added 
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find leads to very precise determinations of /3's; in fact, the 1 — o fitting uncertainties 
are quite modest, 5(3 ~ 0.05 — 0.03/?. Let us notice that the present method does not 
assume a "maximum disk" solution, rather, it is intimely related with the idea that the 
mass distribution in spirals follows the regime of Inner Baryon Dominance proposed by 
Salucci & Persic, (1999b) and then supported by Ratnam & Salucci, 2000 and Salucci et al, 
2000. Finally, the present method computes independently of the value of the disk mass- 
to-light ratio whose uncertainty, [§], large in DM-dominated objects, is therefore irrelevant 
for the aim of this work. 

The disk mass fraction (3 correlates tightly with V (see Fig (2)). In addition, at a fixed 
V opt , part of its scatter originates from observational errors and then it is unrelated to a 
cosmic variance of the halo mass distribution. In any case, conservatively, we derive the 
r.m.s. of the above relation by performing the usual least squares fit 

= 0.75 - 0.95V r.m.s. = 0.05 (9) 

A relationship like eq. (9) has been found from the mass modelling of (smaller samples of) 
spirals (e.g. Persic and Salucci 1988; 1990a; 1990b); it can be considered as the basic law of 
the dark-luminous coupling in galaxies (see Salucci, 1997). 

We derive the values of for the objects of our sample by setting V; um = —0.24 
and inserting in eq (8) the corresponding values of (3 and V. In Fig (4) they are plotted 
as a function of logV op t'. we immediately realize that V/i is roughly constant over the whole 
sample and it shows no systematic variations along the luminosity sequence. Such variations, 
if present, should have clearly appeared given the high-precision measurement of In fact 
the uncertainties 5(3 < 0.04 on (3 propagate into eq (8) in a modest way, as it can be realized 
by differentiating eq (8) and combining it with eq (9) to get 5Vh — (0.25 + V)^ 1 ^ < 0.1 

The average value found for < >= 0.8 ± 0.06 indicates that, at R op t, the halo RC 
is steeply increasing, marginally compatible with a solid body law, Vh oc R. Consequently, 
around R op t, the halo density must decrease with radius less steeply than R~ 0A . By assuming 
a pseudo-isothermal density distribution, Ph{R) (R 2 + a2 ) _1 , it follows that a > 1.3R opt , 
i.e. a core radius significantly larger than R^, the size of region where the bulk of the stellar 
component is located. 

We stress the robustness of this result by noticing that, from eq (8), crucially lower values 
of Vh (i.e. Vft < 0.5) are possible only in the particular combination of a "flat" RC, (i.e. 

§ Since M d = 1.1G- 1 /3 d V^ pt R opt this is ~ (J3 + S0)/(fl - 5(3) 
© 0000 RAS, MNRAS 000, 000-000 
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Figure 3. The innermost regions of Universal Rotation Curve of galaxies of different luminosity (points) compared with the 
disk (+bulge) models (solid lines). Details are in Salucci et. al, (2000) 

for V < 0.2) and a "light disk" (i. e. (3 < (0.5 - V)/0.77 = /3 . 5 ). Also this extreme case, 
however, must be ruled out since such a flat RC cannot be fitted by a mass model with 
V/J less than or equal to (/3o.5, 0.5). (Salucci and Persic, 1999b). 

3 DISCUSSION 

We have investigated the mass distribution of DM halos for a large number of spirals in a 
way which is complementary to the mass modeling of DM-dominated rotation curves. This 
method has been applied to RC's of galaxies of all luminosities, including the most luminous 
ones for which the standard mass modeling is quite uncertain. In detail, for each halo we 
have derived a single but most crucial structural property, namely, its velocity slope at R op t- 

The results are impressive: the halo mass profiles at R opt turn out to be i) independent 
of the galaxy properties, ii) Universal and in) essentially featureless in the sense that for 
any spiral the stellar disk is embedded within a constant density sphere. 

High-resolution simulations and analytical studies of the Cold Dark Matter scenario have 
pointed to a universal halo density profile (Navarro, Frenk and White, 1996) Pcdm(R) oc 
R^ 1 (R + R s )~ 2 , with R s being a function of the dark halo mass M 20 o, of the assumed 
cosmological parameters, and of the red-shift of formation of the halo. Moreover, halos with 
identical values of the above quantities, can still have different R s if formed at various red- 
shifts and/or assembled through different merging histories. In detail, the circular velocity 
of a CDM halo is given by: 

Vcdm(x) 2 oc x _1 (-cx/(1 + cx) +ln(l + cx)) 4^ c(Q , M 200 ) < 30 (8) 

© 0000 RAS, MNRAS 000, 000-000 
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Figure 4. The dark halo slopes Vj, as a function of V op t- As a comparison, in CDM —0.1 < V/j < 0.5 



where c, is the concentration parameter, x = R/R200, and i?2oo(^2oo, z ) is the halo virial 
radius defined by 4/37r.Rij 00 2 00p l A)(l + zf = M 200 - 

Before proceeding further let us notice that we will test the CDM halos at R op t, where 
the baryon infall has not significantly altered the original DM halo velocity profile: R op t is 
external to the region into which most of the baryons have collapsed (Blumenthal et al, 
1986). Violent dark halos-baryonic matter couplings, such as those proposed by Navarro 
et al, (1996) and Gelato and Sommer-Larsen, (1998) can instead modify the original halo 
distribution everywhere; however, given the very heuristic nature of these processes, it is 
best to first compare the standard- infall CDM halos with the galactic halo, and then to 
consider the possibly emerging discrepancies in terms of new theoretical scenarios. 

The highest possible value for ^/^ DM is 0.5, that is achieved on the ~ lOkpc scale only 
for c < 5 (see Bullock et al. 1999 and Navarro, 1998), i.e. for low values of the concentration 
parameter, a property of low-Q universes. This value is quite inconsistent with the average 
value found in spiral dark halos, especially if one considers that high resolution N-body 
simulations converge to a maximum value of \/^ DM = 1/4 (Moore et al, 1998). 

Of crucial importance is also the absence of a significant scatter in the vs. logV op t 
relationship. In fact, the CDM theory predicts that, in a very wide region centered at ~ lOkpc 
and certainly including R opt independently of its relation with the virial radius, galactic halos 
with the same mass do not follow a unique velocity curve but a family of them. These can 
be described by a set of straight-lines with slopes varying between —0.1 and +0.5 (e.g. see 

© 0000 RAS, MNRAS 000, 000-000 
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fig 6 of Bullock et al, 1999). According to CDM the — logV opt plane should be filled 
well beyond the tiny strip of Fig (4). Taken at its face value, the observational constraint 
variance( < 0.1) could imply, within the CDM scenario, that protospiral halos are co-eval 
and have similar merging histories. A second possibility may be that the disk lenght-scale 
R op t, in units of virial radius, is strongly coupled with the structure of the dark matter halo, 
(e.g., Mo, Mao & White 1998; Dalcanton, Spergel & Summers 1997; van den Bosch et al. 
1999, but see also Bullock et al, 1999). 

The sizes of the DM core radii turn out to be very large, at least (3 — 4)i?^, and in- 
dependent of galaxy luminosity and DM mass fraction. This may pose a problem for the 
suggestion that they were formed through some luminous-dark dynamical coupling and it 
may call or for a primordial origin for these "warm" regions of constant density embedding 
the luminous matter or for some dissipative self-regulating" process. 

Burkert (1995) and Salucci & Burkert (2000), from the analysis of individual RC's and 
of the URC, have proposed that the DM halos around disk systems follow the distribution: 
Pb(R) oc (R + Rb)~ 1 (R~ 2 + R~b 2 )i i- e - an NFW profile at large radii which converges to a 
constant value at inner radii. The core radius, Rb is found to increase from 5 kpc to 30 kpc, 
as V opt increases form 75km/s to 300A;m/s. This implies f or the DM halo: — 0.73 almost 
independently of V opt and in very good agreement with the halo slope determinations of fig 
(3). 

Then, since the derived density distribution of galaxy halos is quite different from the 
theoretical one out to the outermost velocity data available (see also Corbelli and Salucci, 
1999; Salucci & Burkert, 2000), we should seriously begin to consider the possibility that 
cosmological processes have cut the link between the initial conditions and the present-day 
galaxies properties. 
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